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Abstract. 

We analyze the form of the probability distribution function P„ (w) of the 



of a given n x n /3-Gaussian random matrix, /? being the Dyson symmetry index. 
This variable, by definition, can be considered as a measure of how any individual 
eigenvalue deviates from the arithmetic mean value of all eigenvalues of a given random 
matrix, and its distribution is calculated with respect to the ensemble of such (3- 
Gaussian random matrices. We show that in the asymptotic limit n — > oo and for 
arbitrary /3 the distribution P„ ,(w) converges to the Marcenko-Pastur form, i.e., 
is defined as P„ (tu) ~ y/(4 — w)/w for w € [0,4] and equals zero outside of the 
support. Furthermore, for Gaussian unitary (/3 = 2) ensembles we present exact 
explicit expressions for P« — (to) which are valid for arbitrary n and analyze their 
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1. Introduction 

Random covariance matrices were introduced by J. Wishart in his studies of multivariate 
populations p] . In physical literature, statistical properties of the eigenvalues of random 
matrices have attracted a great deal of attention since the seminal works of Wigner [2] , 
Dyson [3j H] and Mehta [5]. Various random variables associated with eigenvalues of 
random matrices have been analyzed, such as, e.g., gaps in the eigenvalue spectra, 
number of eigenvalues in a given interval, largest or smallest eigenvalues and etc, with 
a special emphasis put on their typical or atypical behavior. A variety of results and 
their relevance to physical systems have been recently discussed in Ref . [6j [7] . 

One of such variables is the so-called Schmidt eigenvalue, used to characterize, e.g., 
the degree of entanglement of random pure states in bipartite quantum systems. It is 
defined as one of the eigenvalues of a given random matrix divided by the trace, i.e., the 
sum of all eigenvalues. On physical grounds, this variable can be therefore considered as 
a measure of heterogeneity of the eigenvalues and shows how any individual eigenvalue 
deviates from the arithmetic mean of all eigenvalues of a given random matrix. A number 
of significant results on the distributions of such eigenvalues and their extreme values 
for /3-Laguerre- Wishart matrices have been obtained (see, e.g., Refs. [8| |9| [TO j [TT j IT2] and 
references therein). Such random variables have also been considered recently within a 
different context as probes of an effective broadness of the first passage time distributions 
in bounded domains [131 E 



In this paper we analyze the forms of the probability distribution function Pn (w) 
of a Schmidt-like random variable 

Ei=l x v n 

where Xj are the eigenvalues of a given n x n /3- Gaussian random matrix and (3 is 
the Dyson symmetry index. Note that we use a term " Schmidt- like random variable" 
since here we define w as the ratio of a squared eigenvalue over the sum of all squared 
eigenvalues. Within such a definition w is always positive definite and has a support on 

[o,4 

The probability distribution function Pn\w) is given by 

*»M = (*(»-f|)), (2) 
where the average is to be calculated with the weight 




P(xi,x 2 ,...,x„) = -r^exp -- \ x k \ I \xj - Xtf , (3) 



j>i 



K n being a known normalization constant [5]. Our aim is to determine an asymptotic 
behavior of Pn l3 \w) for arbitrary (3 and n — > oo. Apart of this, we will present an exact, 
explicit results for Gaussian Unitary ((3 = 2) Ensembles (GUE) valid for arbitrary n. 

The paper is outlined as follows: In section [2] we provide some general results for 
Pn\w) and analyze its asymptotic forms when n — > oo. In section [3] we present explicit 
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results for the distribution function of the Schmidt-like eigenvalues for GUE. Finally, in 
section H] we conclude with a brief recapitulation of our results. 

2. Asymptotic behavior for arbitrary (3 

Taking advantage of the Fourier cosine representation of the delta-function, we can 
conveniently rewrite Eq. (J2J) as 

P^W = IJ^dycasiwy) (cos (gf)) , (4) 
so that, expanding the cosine into the Taylor series, we obtain 

where pn(x) is the eigenvalue density. Next, using the integral identity 

\ —2k i r°° / \ 



r(2fc) 

we can straightforwardly calculate the multiple integrals over dxi with % = 2,3, ...n, 
which yields 



where fp = n + j3n(n — l)/2. Further on, performing the summation over k in Eq. (JTj), 
we obtain 

P w {w) = YiiiB. r dy COS ( Wy) r dx p (f) (x) {i u/2 _ x (x^j + ^ ( 8 ) 

t/ «/ — oo 

where = (.z/2) ""I,, (2), being the modified Bessel function, and "cc" stands 

for the complex conjugate of l u (z). Equation constitutes our main general result 
valid for arbitrary (3. 

The result in Eq. (jSJ) allows us to establish, for arbitrary /3, the limiting asymptotic 
behavior of the distribution when n — > oo. To do this, we first replace the function I v (z) 
by its integral representation [15] 



and notice that for large values of n the exponential part in the integrand in Eq. (Q is a 
strongly oscillating function of the argument z. This permits us to make the replacement 
(1 — t 2 )^- 1 / 2 pa exp (— z/t 2 ), such that after the substitution w = £a/z/, Eq. becomes 
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Substituting the latter equation into Eq. (JSJ), we arrive at the following 
representation 

Pn\™) = ' r ^/ 2) 1W01 Hdycosiwy) f°° dx (x) cos , (H) 

7r v//8/2r[(/^ - l)/2] Jo J-oo \ IP J 

which yields, after performing the integrations, the following asymptotic form 



i.e., it simply expresses the desired probability distribution of the Schmidt-like random 
variable w through the eigenvalue density with an appropriately rescaled variable. The 
asymptotic behavior of the latter is well-known and is defined by the Wigner semi- 
circle distribution [21 [5] , so that after some very straightforward calculations we find the 
following asymptotic form for the normalized probability distribution function : 



27T { 0, for w > 4 . 

Equation (1131) holds for any value of the Dyson symmetry index /3. It might seem 
surprising at the first glance that the limiting distribution in Eq. ffTB"]) has the form of 
the Marcenko-Pastur law [IB]. On the other hand, recall that asn-> oo, the eigenvalues 
tend to be equidistantly-spaced so that the sum ^2jX 2 /n tends to a constant. Then, it 
becomes clear why the distribution Pn\w) converges to an appropriately normalized 
single eigenvalue density, defined by the semi-circle distribution [2J [5] , so that its squared 
value is distributed according to the Marcenko-Pastur law. 

3. Gaussian Unitary Ensemble 

We turn now to the GUE case (/3 = 2 and fa = n 2 ), aiming to evaluate an explicit 

(2) 

expression for the probability distribution Pn (w), valid for an arbitrary value of n. In 
this case, the eigenvalue density is given by 

= ^fr? - H n+1 (x)H n ^(x)} , (14) 



where H n (x) denotes the Hermite polynomial [15]. Using Eqs. (Q and (I14p . we can 
represent the integral as 

/oo rl ^ n _ ^2-j( n2 - 3 )/ 2 e inyt 2 

j X pV\x)i h/ ^ {2x y /i^i)=J_ i 2nnlr(n2/2 _ 1/2)7r HtV^y) , (is) 

where 

/oo 
dxe - {x - tV ^? ^2 (x) _ Hn+1 ( x)Hn ^ x) } . (i 6 ) 
-oo 
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Further on, this function can be expressed in terms of the associated Laguerre 
polynomials (p) since for n > m [T5j 

/oo 
dxe- (x - z)2 H m (x)H n (x) = 2 n ^m\z n - m Ll n - m \-2z 2 ) , (17) 
-oo 

which leads to 

f(ty/i^) = 2" v ^4 1 2 1 (-2t 2 zn 2/ ), (18) 



where we made use of the following recurrence relation between the associated Laguerre 
polynomials : nL n (p)+pL^}_ 1 (p) = L^_ x (p). Then, the probability distribution function 
Pn (w) becomes 

dycos(wy) J ^= (1 - v ) (n2 - 3)/2 (e inyv L^^inyv) + ccj , (19) 
where the normalization constant C is given explicitly by 

nr[(n 2 - l)/2] ' 1 ' 

and "cc" here stands for the complex conjugate of exp (inyv) L^ l 1 2 1 (— 2inyv). Using next 
the series representation of the Laguerre polynomials, the integration over dy reduces 
to the calculation of the following integrals 

n— 1 /»cxd 

VflJ dycos(wy) e inyv (-2viny) k + e- inyt \2viny) k , (21) 



or 

n-l 



2 ]P ak J d y cos ( w y) cos (ny v + ~^-J (-2vny) k , (22) 

which, using the fact that cos(a; + kit/ 2) = cosx, can be put into the following form 
nL^J^v-^jSiw-nv). (23) 

Next, the integration over dv can be performed by parts, taking into account that 
yfc-i/2^j _ t ,yn 2 -3)/2 with fa > q van ishes at the integration limits. Then, the operator 
expression 



r(n 2 /2) r(1) Id \ (1 - ■ 1 ,)<" 2 - 3 " 2 



0Fn 2 r[(n 2 -l)/2] V / v 7 ^ 



(24) 

v=w/n 



is obtained where powers of the polynomial operator are understood as 
Explicitly, the action of the Laguerre polynomial operator is defined as 

L m ( 2 d\ g-vr-^ f ( -2)>ni •LL^nn S"-*>") (25) 
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where the derivatives can be identified with the Rodrigues formula for the Jacobi 
polynomials Pn ,b \p) P3] as 

J_ (>-l/ 2 (1 _ v) (n>sy^ = (~l) fc fc! (1 -v)^-W- k p ^-S )/2 - k ,- im{2v _ 1} (26) 

Consequently, recalling that v = w/n, we find the following explicit result for the 
distribution P n (w) : 

P m M = r( B y2) Mj^ x 



x 

fc=0 



2*(l-^""-pr-'™>(2^-l). (27) 



Further on, the sum on the right-hand-side of the latter equation can be also represented, 
after some straightforward calculations, as a polynomial of w, which yields 

P (2) H _ 2r(n 2 /2) (l-w/nt 2 - 2 ^ 2 ^ 



7r3/4r((n 2 - l)/2)nV2 ^ 

x£<=£(»-V»" (28) 

m=0 v 



where the coefficients defined by 

a m = f (1 - 2t)— 1 2 F, f-m, ^ - 1,1 2t) , (29) 

2-Fx (•) being the Gauss hypergeometric function. Equations (I2"7|) and (I2"81 constitute 
our principal results for case of the Gaussian Unitary Ensemble. 

Before we turn to the analysis of the asymptotic behavior of the distribution 
function, it might be expedient to present several first Pn (u>) explicitly. Below we 
display Pn\w) for n = 2, 3, 4, 5 and 6: 

P?\w) = \^=L=, (30) 

7T y/{2 -W)W 



P?\w) 



35 (3 - w) 
576^3 



(3-2w + 3w 2 ), (31) 



P 4 (2) M = ^—^L ( 12 + 30 w - 53 w 2 + 38 w 3 ) , (32) 
1716 7r 



and 



2028117 5- w 7 / „ „ a , , 

i ^— 375 - 300 w + 4490 w 2 - 5996 w 3 + 2711 w 4 , 33 

81920000000000 y/E y/w V / 
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32768 (6 - w)™/ 2 

' 810 + 3780 w- 18090 w 2 



25113523969051155 irJw 



52878 w 3 - 49567 w 4 + 16144 w 5 



(34) 



One notices that the expressions in Eqs. ( 1301 to (134j) all diverge as 1/ a/w when 
u; — ?• 0. Next, all these expressions vanish (for n > 2) as a power-law at the other edge 
of the support [0, n], with an exponent dependent on n. The case n = 2 is special: 
P^. 2 ( w ) diverges at both edges and has a minimum at w — 1, which signifies that in 
2x2 Gaussian random matrices the two eigenvalues are most probably very different 
from each other. Further on, in Fig. ([T]) we plot these explicit forms together with 
more lengthy expressions for n = 7 and n = 12. One observes that for w < 1 the 
distributions of arbitrary order are very close to the asymptotic result in Eq. ( fTBT) . The 
distibutions are multimodal indicating a set of probable and unprobable values of w, 
which mirrors certain structuring of the eigenvalues. As n gets progressively larger, 
the distributions become closer to the asymptotic result, Eq. ( TT3"|) . for any w G [0,4]. 
Curiously enough, despite a rather complicated form of the polynomials in the second 
line of Eqs. f|2T|) and (|2"8|) . they all show an appreciable variation with w only for w < 4 
and are indistinguishable from zero for larger values of w, despite the fact that formally 
their support extends to larger than 4 values of w. 




w 



( 2) 

Figure 1. (color online) The distribution P„ (w) in Eqs. (|27l) and ((28)) for n = 2 
(blue), n = 3 (green), n = 4 (purple), n = 5 (pink), n = 6 (magenta), n — 7 (orange) 
and n = 12 (red). Solid black line defines the asymptotic result in Eq. (|13j) . 



For arbitrary n, an asymptotic behavior of Pn 
can be readily deduced from Eq. (|2 



(2), 



w) for m; 1 and w close to n 



As we have already remarked, one finds that for 
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w 0, the distribution shows a generic singular behavior of the form: 

^. M^ ', (35) 
y/wnT ((n 2 - l)/2) 

where the amplitude 

r(n 2 /2) ^(-71+1,1,2, 2) ^ 1 



0mr((n 2 - l)/2) 7T 

when n — > oo, in agreement with the general result in Eq. ( fl~3|) . This implies, in turn, 
that fo a given random matrix a randomly chosen eigenvalue will most probably be much 
less than the arithmetic mean of all eigenvalues. Further on, on the opposite extremity 
of the support, when w is close to n, we have from Eq. ( |2"7|) that 



n(2)/ \ 2™ 1 T(n 2 /2) x(n 2 -2n-l)/2 / or .\ 

P - H ~ ^(n-l)\ n(" 3 -^r((n-l)V2) ( " " " } ' ^ 

i.e., Pr(w) attains a zero value as a power-law when w — >■ n with an exponent which 
grows in proportion to n 2 when n — >• oo. This implies, in turn, that for w sufficiently 

(2) 

close to n the value of Pn (w) decays faster than exponentially with n. 

Finally, we address the question how the Marcenko-Pastur law in Eq. (fl~3|) can 
be derived from our Eq. (1271) . Below we briefly outline the steps involved in such a 
derivation. Note first that for w < n, one has 

W\ (" 2 -2n.-l)/2 / nW \ 



n 

as n — y oo, so that 



^exp(-— j, (38) 



■ 

> 



T(n 2 /2) (l-w/n) (n - 2 "" 1 )/ 2 exp (-nw/2) 

?r((n 2 -l)/2)n3/2 ^ ^ ' ( } 



Further on, one has that, as n — > oo, [15 



^^-^o-iWv?)- < 4o) 

Using next the integral representation of the Hermite polynomials 



2 / y7T J Jo V 

we can resummate the series in Eq. ( |2"7|) to find that, as n — > oo, the sum in the latter 
equation converges to 

n-l 



fe=0 v 7 



Note next that as n — > oo US 



^ ( l_\ ^ n^J, (2y/Vg ^ (43) 
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so that the integral on the right-hand-side of Eq. (|42p converges to 

Vr-»(r)f *«*(-£) w 

where J\(x) is the Bessel function. Noticing finally that 

^ (-2^) - 1 < 45) 
as n — > oo, we can perform the integral in Eq. ( jUJ) in this limit, to get 

exp « { V 71 " for < ^ < 4, 

V vr P V 2 / 1 0, for w > 4. V ; 

On combining the latter equation with Eq. (139|) . we arrive at the Marcenko-Pastur law 
in Eq. . 



4. Conclusions 

To recap, we analyzed the probability distribution function P„ («;) of the Schmidt-like 
random variable w = x\ / \Y^ =1 x j/ n J j Eq. ([ID, where are the eigenvalues of a given 
n x n /3-Gaussian random matrix. This variable, by definition, can be considered as a 
measure of how any individual eigenvalue deviates from the arithmetic mean value of all 
eigenvalues of a given random matrix, and its distribution is calculated with respect to 
the ensemble of such /3-Gaussian random matrices. We showed that for arbitrary Dyson 
symmetry index /3 in the asymptotic limit n — > oo the distribution Pn\w) converges 
to the Marcenko-Pastur form, i.e., is defined as PP(w) ~ y/(4-w)/w for w G [0,4] 
and equals zero outside of the support. For Gaussian unitary (j3 = 2) ensembles we 
presented exact explicit expressions for Pn l3 ~ 2 \w) valid for arbitrary n. We realised 
that, in general, Pn~ 2 \w) has a multimodal form indicating probable and unprovable 
values of w, which mirrors certain structuring of the eigenvalues. We realised that the 
convergence to the Marcenko-Pastur form is rather fast, so that already for n = 12 the 
exact result appears to be quite close to the asymptotic form. 



Acknowledgments 

The authors acknowledge helpful discussions with Oriol Bohigas and Satya N Majumdar. 
This work is supported by the Brazilian agencies CNPq and FAPESP. GO is partially 
supported by the ESF Research Network "Exploring the Physics of Small Devices". 



References 

[1] Wishart J (1928) Biomctrika 20 A 32 
[2] Wigncr E P (1951) Proc. Cambridge Philos. Soc. 47 790 
[3] Dyson F J (1962) J. Math. Phys. 3 140; 3 157; 3 166 
[4] Dyson F J and Mchta M L (1962) J. Math. Phys. 3 701 



9 



Distribution of Schmidt-like eigenvalues 



M. P. Pato and G. Oshanin 



[5] Mehta M L (2004) Random Matrices, (Academic Press, 3d Edition, London). 

[6] Bohigas O and Pato M P (2010) J. Phys. A 43 365001 

[7] Majumdar S N, Nadal C, Scadicchio A and Vivo P (2011) Phys. Rev. E 83 041105 

[8] Lloyd S and Pagcls H (1988) Ann. Phys. 188 186 

[9] Kubotani H, Adachi S and Toda M (2008) Phys. Rev. Lett. 100 240501 

[10] Akemann G and Vivo P (2011) J. Stat. Mech. P05020 

[11] Vivo P (2011) J. Stat. Mech. P01022 

[12] Mejia-Monasterio C, Oshanin G and Schehr G (2011) Phys. Rev. E 84 035203 

[13] Mejia-Monasterio C, Oshanin G and Schehr G (2011) J. Stat. Mech. P06022 

[14] Mattos T G, Mejia-Monasterio C, Metzler R and Oshanin G (2012) Phys. Rev. E 86 031143 

[15] Abramowitz M and Stegun I (1965) Handbook of Mathematical Functions (New York, Dover). 

[16] Marcenko V A and Pastur L A (1967) Math. USSR-Sb. 1 457 



10 



